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Abstract. We derive an upper bound on the action of a direct product of two 
quantum maps (channels) acting on multi-partite quantum states. We assume that the 
individual channels Aj affect single-particle states so, that for an arbitrary input pj, 
the distance Dj{Aj[pj], pj) between the input pj and the output Aj[pj] of the channel 
is less than e. Given this assumption we show that for an arbitrary separable two- 
partite state pi2 the distance between the input pi2 and the output Ai ® A2[pi2] fulfills 

the bound i:)i2(Ai (g) A2[pi2],pi2) < ^^2 + 2^(1- l/di)(l - l/o?2) e where di and 
are dimensions of first and second quantum system respectively. On the contrary, 
entangled states are transformed in such a way, that the bound on the action of the 
local channels is -Di2(Ai (8) A2[pi2], P12) < 2-^2 — l/d e, where d is the dimension of 
the smaller of the two quantum systems passing through the channels. Our results 
show that the fundamental distinction between the set of separable and the set of 
entangled states results into two different bounds which in turn can be exploited for a 
discrimination between the two sets of states. We generalize our results to multi-partite 
channels. 
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1. Introduction 

Investigation of properties of communication channels is more than ever today a central 
issue of information sciences. It is generally accepted that quantum systems have the 
capacity to carry information efficiently and any transformation of these systems can 
be considered as an action of a quantum channel (see e.g. Refs. [H El El HI [5] ) . 

Some general questions arising from the transmission of quantum entaglement 
through quantum channels has been analyzed by Schumacher in Ref. [6]. He has 
considered a pure entangled state of a pair of two systems R and Q and the system 
Q has been subjected to a dynamical evolution (quantum channel). Schumacher has 
shown that the two quantities of interest, the entanglement fidelity and the entropy 
exchange Se, can be related to various other fidelities and entropies and are connected 
by an inequality reminiscent of the Fano inequality of classical information theory. 

In this paper we address the question how two local channels (each acting 
independently) do affect a bi-partite quantum state. This scenario is rather general 
and can be applied to a number of situations, e.g. quantum computation with quantum 
computer imperfectly isolated from environment or analysis of quantum error correcting 
codes. In the context of quantum error correction this problem has been addressed by 
Knill and Laflamme in Ref. [7] for a particular case of two qubits and for a particular 
choice of a distance (fidelity) characterizing the change of the bi-partite state. In Ref. [8] 
Aharonov et al. have analyzed errors for a general model of quantum computation 
with mixed states and non-unitary operations. There however a different measure was 
introduced. The rationale being that measurable distinguishability of gates (super- 
operators) should not increase if we consider additional quantum systems which do 
not evolve. Here on the contrary we are not interested in the distinguishability of 
superoperators but rather in the actions of the channels on a given state and how to 
relate these local actions with the change of the global state. 

Specifically, consider a pair of quantum channels characterized by maps Ai and 
A2, respectively. It means that after sending a quantum system over, for instance, the 
first channel the final state of the quantum system (or equivalently the output of the 
channel) is Ai[pi] where pi is the corresponding input. Moreover, let the two channels 
fulfill the following condition 

D,(A,[p,], p,) < e, \/pj G Sin,), J = 1, 2 , (1) 

where Dj{.,.) for j = 1,2 are some distance functions (metric) defined on the set of 
all density operators S(Hi) and 5 (7^2) representing the set of all physically realizable 
states of quantum systems passing through the channels 1 and 2, respectively. These 
conditions restrict the action of each of the two channels independently of the action 
of the other channel. Specifically, the state of a quantum system affected by one of the 
two channels has to be in a small (epsilon) neighborhood of the state describing the 
quantum system before the system was sent through the channel. 

The parameter e quantifies the action of the two quantum channels. For e = 
the two channels are "perfect" (i.e., noiseless, that is, the information transmitted via 
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channels is not disturbed) as the output equals to the input while for e large the output 
can be significantly different from the corresponding input 0. 

The question we would like to address is, how big the change induced by the two 
local channels is when the inputs are correlated. That is let us prepare an arbitrary 
initial state pi2 of two quantum systems. The first part of the jointly prepared system 
is sent over the first channel while the second part is sent over the second channel. Both 
channels individually fulfill the condition ([1]), where, e.g. pi = Tr2Pi2- In what follows 
we will show that the two-partite action of the channel Ai ® A2 for all possible physical 
states P12 G iS(7Yi2) fulfills a bound on its action that is determined by single-partite 
conditions given by Eq. ([T]). 

Let us note that the problem can be transformed into the estimation of the map 
= Ai ® A2 — I12 where the map I12 is the identity acting on the joint system. If the 
distance -Di2(-, •) as well as distances Di{., .) and D2{., .) are defined via a norm then 
our task is to estimate the norm ||f2(pi2)||. Similar expressions for a general class of 
the so called p-norms has been studied extensively for Q being a physical map (more 
specifically the product of two physical maps) in Refs. [9], [10]. However, in our case the 
map Q is neither a positive map nor a direct product of two maps. Due to the fact that 
the map fl is not positive and subsequently not physical our situation is not applicable 
to Refs. [9l [To] and similar studies. 

The paper is organized as follows. In Sec. 2 we introduce necessary definitions and 
discuss a particular case of separable states, i.e. the initial state of the joint system 
(the system composed of two quantum systems that are sent over the two quantum 
channels) is separable. As a next step we drop any assumptions on the initial state and 
analyze the most general case of an arbitrary initial state in Sec. 3. The results obtained 
are discussed in Subsection 3.1. In Sec. 4 we extend our analysis to the case of more 
than two quantum channels and illustrate the nature of changes on a simple example. 
Finally, in Sec. 5 we summarize our results and outline possible extensions. 

2. Separable inputs 

In the formulation of the problem we encounter three different metric (distance) 
functions: Di{.,.), D2{.,.) and Di2{.,.) acting on different sets and thus measuring 
distances between different types of objects. In order to make our discussion explicit we 
will consider a specific choice of the distances offered by the norm of the Hilbert-Schmidt 

I Let us note that there is no relation between the parameter e and the capacity of the channel in 
general. 
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spaces corresponding to the systems 1, 2 and the joint system 12, respectively [§| 

DaiPa,Cra) = \ \pa - (^a\\a 

= ^jTTa[{pa-aa){pa-CTay] • (2) 

The label a denotes the system 1, 2 or the joint system 12 and pa, Ca G SiTia) are density 
operators representing possible physical states of the system labelled a. The norms that 
we have used to define distances _Di(.,.), _D2(.,.) and Di2{.,.) are called 2-norms and 
are only a particular case of the so called p-norms. However, due to the fact that we 
will use only basic properties of the distances .), D2{-, ■) and -Di2(-, •) we will keep 
our derivation as general as possible so that it can be repeated with a broad class of 
different distances. Only in the end we will use the specific choice of distances to derive 
a tight bound. 

Our task is to estimate the distance 

Di2(Al ® A2[P12],P12) , (3) 

for all physically reasonable initial states pi2 G 5(7ii2) provided the two maps Ai and 
A2 fulfill the condition ([T]). First, note that for any distance inequality (this follows from 
the triangle property of a distance) holds 

Di2(Ai ® A2[pi2], P12) < ^12(Al ® l[pi2], P12) 

+ Di2(l ® A2[Ai ® l[pi2]].Ai ® l[pi2]) . (4) 

It means that instead of considering the case with two local channels it is sufficient to 
consider only an action of a single local channel acting on one of the two subsystems 
and estimate the distance Di2(Ai ® l[pi2])Pi2)- 

We start with the simplest case - the case of factorizable states of the form 
P12 = Pi® P2- This corresponds to the situation as if the two channels were considered 
separately so that the two quantum systems that are sent through the channels are 
prepared individually. In this case we exploit the following property 

^i2(pi ®P2,Pl ® P2) < Di{p[,pi) , (5) 

of the distances Di{., .), D2{., .) and D12 where the operators pi, p[ and p2 are density 
operators representing states of the first and the second system, respectively. Let us 
note that this relation holds even if the distances are defined with any p-norm or even 
fidelity. Using Eq. ([5]) we have that Di2(Ai ® l[pi ® P2]5Pi ® P2) < -Di(Ai[pi], pi) and 
consequently, for the initial state of the form pi® P2, the distance ([3]) is always less or 
at most equal to 2e 

Di2(Ai®A2[pi®p2],Pi®P2) <2e, (6) 

§ The set of all density operators representing the set of physical states of a quantum system is a 
subset of a vector space. In such case it is natural to define the metric (distance function) with the 
help of a norm so that the linear structure of the vector space is respected. There are several ways 
how to introduce a norm on a vector space. However, the set of all density operators is also a subset of 
the Hilbert-Schmidt space which is a Hilbert space and we can use the norm induced with the scalar 
product of the Hilbert space. 
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due to Eqs. (Hj) and ([I]). 

The same holds for the initial state pi2 of the form pi2 = J2i (^iPi®p2i where > 0, 
«j = 1 and p\ and denote density operators of the system 1 and 2 respectively, 
that follows from the linearity of the map Ai ® A2 

Di2(Ai ® A2 E "iPi ® E "iPi ® P2) 

i i 

j i 

and the fact that the distance -Di2(-, •) is jointly convex , that is 

Di2(E "j^i ® ^2 [pi ® pi], E "iPi ® P2) 

j j 

< E ai^i2(Ai ® A2[pi ® pi], pi ® pi) . 

i 

The last expression is a sum of terms where each term is bounded by 2e and the sum of 
the coefficients is equal to unity. In consequence we obtain the bound 

Di2(Ai (g) A2E "iPi ® P2]. E "»Pi ® P2) < 2e , (7) 

i i 

for an arbitrary separable state. 
2.1. Hilbert- Schmidt distance 

The bound on the action of a product of two quantum channels on separable states ([7]) 
is vahd for any triple of distances Di{., .), D2{., .) and -Di2(., .) that satisfy relation 
(the distance D2{.,.) has to fulfill the relation ([5]) with swapped labels 1 and 2) and 
in addition the distance -Di2(-, •) has to be jointly convex. That is, the bound is valid 
if D2{.,.) and Di2{.,.) are trace distances [Jj] or, more generally, the distances 

defined with p-norms or even fidelity. The question is whether it is possible to derive a 
better (tighter) bound or, in other words, whether the bound is optimal. For the trace 
distances the bound is optimal indeed and it can be shown that there is a pair of maps 
such that the bound is saturated. In what follows we will show that for the distances 
introduced in Eq. ([2]) the bound can be further optimized. 

Let pi = 1/dil+c.a be an input of the channel 1. We have expressed the state of the 
system labeled as 1 using the identity operator 1 and df — 1 generators a = {ai, 0-2, • • •} 
of the group SU{di) multiplied with the complex unity where di is the dimension of 
the Hilbert space of the system 1 and the vector c = {ci, . . .} is a real vector with 
d1 — 1 elements. In addition we require that the set of operators {era} satisfy the 
ortho-normalization condition Tr aaCTf^ = dap. After the quantum system has been 
sent through the quantum channel Ai the state of the system (the output) can be 
expressed using the same notation Ai[pi] = l/dl + c' .a with new coefficients c' where 
the prime indicates the fact that the state has been sent through the quantum channel. 

II The trace distance is defined with the help of the l-norm and D{p, a) is equal to the sum of eigenvalues 
of the positive operator \p — a\ where |p — cr| = ■\/(yO — cr)^(p — u) and p,<J ^ B{H). 
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Equivalently, p2 = l/d2l + d-f is the most general state of the system 2 where f are 
generators of SU{d2) multiphed with complex unity, d2 is the dimension of the Hilbert 
space of the system 2 and the operators {r^} satisfy relation Tr r^gr^ = 5^^}- 

We estimate the distance ([3]) for an arbitrary separable state and the particular 
choice of distances ([2]). Due to the joint convexity of the distance Di2{.,.) and the 
linearity of the map Ai ® A2 it is sufficient to consider the case where the state pi2 is a 
pure state (for more details see the end of the previous section) 

P12 = (l/rfil + c.a) (g) (l/dsl + ^.r) (8) 

where c.c = (1 — 1/di) and hi = (1 — 1/^2)- 

In this case we do not use the relation Eq. (jlj) which means that the two channels 
are not considered separately and the output of the product of the two channels Ai and 
A2 is 

Ai ® A2[pi2] = + c.a) (l/rf2l + b'.f) . (9) 

Inserting the two expressions, input ([8]) and output ([9]), into the definition of the 
distance ([2]) we obtain that 

L)i2(Ai (g) A2[pi2],Pl2) = 

I |(c' - c).a (g) 1/^2! + 1/dil (g (b' - b).f + c'.a (g b'.f - c.a (g fe.f | |i2 . (10) 
Last expression squared can be bounded from above by a sum of three terms 

||(c'-c).a® l/d2l||i2 + ||1M1® (&'-6).r||^2 
+ [ ||(c' - c).a 0b'.f\\i2 + \\c.a0 (b' -b).f\\i2 ]^ • 

Observing that ||(c — c').a (g I/C/2III12 = 1/^2 -Di(Ai[pi], pi), and equivalently ||l/ciil(g 
(b — 6')-^lli2 = l/'^i-D2(A2[p2], P2) and b'.b' < (1 — l/c?2) we can bound the distance 
squared with the expression l/d2e'^ + 1/die^ + (y^l — 1/c/i + y'l — 1/(^2)^6^. Finally, the 
distance between the input and the corresponding output of the product of the two 
channels fulfills the bound 

Di2(Ai®A2[pi2],Pi2) < y2 + 2^(l-l/rfi)(l-l/d2) e, (11) 

where di and ^2 are the dimensions of the Hilbert spaces corresponding to the quantum 
systems sent through the channels 1 and 2, respectively. Even though we have proved 
the bound for pure separable states, we note that the result is valid for an arbitrary 
separable state due to the linearity of the map Ai (g A2 and the joint convexity of the 
distance -Di2(-, •)• The bound (ITT!) is undoubtedly better than the bound ([7]) as it has 
been derived for a specific choice of distances. In addition, it can be shown that the 
bound is optimal in the sense that there is a pair of maps Ai and A2 and a separable 
state P12 such that the bound ( fTTl) is saturated (optimality is discussed in a more detail 
in Sec. |5]). 
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3. Entangled states 

We have seen that if the initial state of the joint system 12 is factorizable or 
even separable then the action of the two channels is bounded by the expression 



2 + 2y(l — l/(ii)(l — 1/^2) e It may be tempting to say that the same holds for an 
arbitrary state. However, as the next example illustrates, if the joint state of the two 
systems 1 and 2 is entangled then for certain maps the separable bound can be broken. 

Let us consider the Hilbert spaces Hi and H2 corresponding to the systems 1 and 
2 to be two-dimensional spaces. This is the simplest possible case though the physical 
representations of such systems are numerous. As an example we can mention spin 
one-half particles, polarized photons or particular internal degrees of freedom of an ion. 
Let us note that in quantum information theory such systems are denoted as qubits 
since they represent quantum analogue of a classical bit of information. 

Then, any physical state of the system 1 (or equivalently of the system 2) can be 
written as pi = ^{l + a.a), where a = (ai, 0:2, as) is a vector in a three dimensional real 
vector space and the three matrices a = (o"i, a2, cr^) are the well known Pauli operators. 
For the matrix pi to represent a physical state the norm of the real vector a has to be 
less or equal to 1. It follows that the set of all physically realizable states of the system 
1 corresponds to a unit ball (Bloch sphere) in the three dimensional vector space IR^. 

The map Ai we will consider in this particular example is a simple contraction of 
the ball representing the set of states such that 

Ai:pi^^{l + il-k)a.a), (12) 

where (1 — A;) is a parameter of the contraction. Physically, the map Ai describes a 
channel with uncolored ( "white" ) noise since each input state is mixed with the absolute 
mixture 1/2 1 which is the fixed point of the Ai. In order to preserve the condition ([T]) 
the parameter k has to fulfill the relation k < ^J2e. In what follows we assume k = V2e. 

In the same way the most general state of the system 2 is p2 = |(1 + P-ff), where 
P = {^1,^2,^3) is a real vector and < 1. The map A2 has been chosen to be the 
same as the map Ai 

A2:p2^^(l + (l-fc')/3.a), (13) 

with the same contraction parameter k' = k = \/2e so that the condition ([T]) is fulfilled 
in this case too. 

To show that the separable bound can be broken we have to consider an entangled 
state. However, we will not consider an arbitrary state but a very specific one - a 
maximally entangled state known as the Bell state of the form pi2 = 1/2(|01) — 
|10))((01| — (10|), where and 1 denote two basis vectors of Hi (or 7^2)- For 
subsequent calculations it is useful to rewrite the state using the Pauli operators 
P12 = 1/4 (1 (8) 1 — o"! CS) cTi — (72 ® ^"2 — ® (Ts). Inserting pi2 into Eq. ([3]) and using the 
linearity of the transformation Ai ® A2 as well as Eq. (I2j) the distance in Eq. ([3]) reads 

£)l2(Ai(g)A2[pi2],Pl2) = II ^ Wl^ai + 0-28)0-2 + 0-3 (8)0-3} 1 1 12 . 
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Both constants, k as well as k' are equal to ^/2e. Neglecting terms of the order and 
evaluating the norm using the scalar product we find 

Di2(Ai ® A2[pi2],Pi2) ~ V6e . (14) 

This result clearly shows that even though the two maps Ai and A2 fulfill the relations ([1]) 
the map Ai ® A2 constructed as a direct product of the two maps can affect the states it 
acts on in a much stronger way. How much the joint (and particularly entangled) states 
can change by two arbitrary maps Ai and A2 is addressed in the next paragraph. 

Let P12 be an arbitrary mixed state. The deviation of the output of the channel Ai (S> A2 
from the input pu is characterized by the distance Eq. ([3]). In order to estimate the 
distance we exploit (as in the case of separable states) the bound given by Eq. (jlj) 

Dl2(Ai®A2[pi2],Pl2] < £'l2(Ai®l[pi2],Pl2) + A2(l®A2[pi2],Pl2) , (15) 

where pi2 = Ai®l[pi2]. As we do not make any assumptions about neither the maps Ai 
and A2 nor the initial state pi2 the two states pi2 and pi2 can be arbitrary physical states 
of the joint quantum system, i.e. arbitrary density operators. It means that taking, for 
instance, the first term on the right-hand side of Eq. (fT5l) we need to estimate this term 
for all possible maps Ai and all possible states pi2. This fact allows us to rewrite the 
bound for ([3]) in a different way 

£)l2(Ai (g) A2[pi2],Pl2) < 2 SUP|A,,pi2}ll^l ® ^[Pl2] - P12II12 , 

where the factor 2 appears because we have two terms in Eq. (ITSl) and the supremum 
runs over all possible maps Ai and all initial states pi2 0. 

A mixed state pi2 can be decomposed into a mixture of pure states pi2 = 
Y^k(^k\i'k){4'k\- Using a basic property of the norm (or joint convexity of the distance) 
and the normalization condition J2k(^k = 1 we can simplify the last expression and 
instead of searching for supremum over all possible states pi2 of the joint system 12 it 
is sufficient to consider pure states only. It means that 

Di2(Ai®A2[pi2],Pl2) < 2 SUP|^^,|^^^^|}|| Ai®l[|^)(7/^|] - II12, (16) 

where the supremum runs over all possible maps Ai and all possible pure states 
€ S(Hi2) of the joint system 12. Since we have used only a basic property 
of the norm the last relation is valid for any distance defined with the help of a norm 
(or more generally any distance that is jointly convex). However, in what follows we 
will use specific properties of the Hilbert-Schmidt norm and further results are valid for 
that particular choice of the norm only. 

Any pure state {ip) G Hu can be expressed using the Schmidt basis 

\ij) = Y.f^k\k)i®\k)2, (17) 

k=l 

% Given the fact that we have not specified the dimension of neither the system 1 nor the system 2, 
the two systems can be different. Therefore we should find supremum over ah Ai's and aU pi2 of the 
first expression in Eq. ([T5|) and aU A2's and aU pi2 of the second expression in Eq. (fT5|) . However, the 
results are the same in both cases. 
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where and {|A;)2} are two sets of orthonormal vectors of Hi and TC2, respectively, 

and jSk are real positive coefficients. The integer denotes the number of elements in 
the Schmidt decomposition of the given pure state and is always less or equal to the 
dimension of the smaller of the two Hilbert spaces Tii and 7^2- In this particular basis 
the state pi2 = has the form 

l^)(^l= E/5/^A|A;)i(/|®|A:)2(/|. (18) 

k,l=l 

Let us now estimate the expression || Ai — II12 from Eq. f|T6|) . 

Using Eq. (ITSl) for the density operator and tracing over the degrees of freedom 
belonging to the second system we have that 

||Ai®l[|V)(^|]-|V^)(^|||L= E PlPfTriVM{V,,)\ (19) 



where = Ai[|A;)i(/|] — \k)i{l\. At this point we apply the relations Eqs. (lA.ip 



(IB. II) and ( IB. 21) (proved in Appendix A and Appendix B ) and Eq. ([T]) that establish the 
following inequalities 

TiiVkiiVkiY <2e^; "ik ^ I ■ 

Tri Vkk{VkkV <e'; VA; . 

These inequalities bound each contribution (trace term) in the sum on the right of 
Eq. (|T9|) . If we replace each term with the corresponding bound and maximize over all 
possible Pj then we do estimate the expression on the left-hand side of the last equality 
as 

II Ai®l[|^)(^|]-|^)(^| ||?2<(2-l/rf)e^ 

where d is the dimension of the smaller of the two Hilbert spaces Tii and 7^2 in case 
the two subsystems 1 and 2 are different Q- Since the result is independent of both the 
map Ai and the state it holds for all maps Ai and all density operators 

(representing pure states). Consequently, the supremum over all maps Ai and all pure 
states P12 is less or equal to this value and so is the distance ([3]) 



Di2(Ai (g) A2bi2], P12) < 2^2 - 1/d e. (20) 
The bound (1201) is valid for entangled as well as separable states. However, for 



separable states we have already found a tighter bound y 2 + 2y (1 — l/di){l — 1/^2)^ 
[see Eq. (fTTl) ] which means that the entangled states can be affected by independent 
channels more strongly than separable states. 

+ If the two subsystems 1 and 2 are different then the number of elements in the Schmidt decomposition 
Eq. (fT7|) is always less or equal to d - the dimension of the smaller of the two Hilbert spaces Hi and 
Ti.2- Consequently, the number of coefficients (3j we maximize over is always bounded by this number, 
which in turn bounds the maximum. 
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3.1. Detection of entanglement 

The difference in the behavior of separable and entangled states resulted into two 
different bounds. The bound for entangled states is weaker and this bound is obeyed by 
entangled as well as separable states. On the other hand the bound for separable states 
Eq. (ITTl) is tighter and need not be fulfilled by entangled states. Subsequently, any state 
that violates the bound ffTTj) is necessarily entangled and a direct product of physical 
channels can be exploited as a kind of "entanglement witness" . Let us point out that the 
entanglement witnesses known in the literature, Refs. pTl [T2] . are based on a different 
approach. They are constructed using positive but not completely positive maps (that 
is non-physical maps) acting on one of the two subsystems and the non-positivity of 
the final operator (output) is the indication of entanglement. On the contrary, in our 
case, we have a product of two physical maps so that a physical (completely positive) 
map is acting on each of the two subsystems and the difference between an input and 
the corresponding output is measured. In addition there is a potential advantage in 
this approach. Not only the question whether a state is entangled or separable can be 
answered. If we relate the distance to entanglement then we could answer the question 
how much entanglement is shared by two quantum systems. 

Similarly as in the case of entanglement witnesses, given a pair of maps, the 
detection need not be (and in general is not) perfect. In other words given a pair 
of channels only a subset of the set of all entangled states violates the bound Eq. ffTTl) 
and those are the only states that are detected as entangled. Naturally, it is desirable to 
optimize the detection so that the whole set of entangled states is detected. There are 
several things we can do to optimize the detection of entangled states using quantum 
channels: 

(i) Optimal choice of the distances Di{., .), D2{., .) and -Di2(., .) 

(ii) Optimal choice of the pair of channels (maps Ai and A2) and subsequent derivation 
of the bound for separable states for that particular choice. 

It is obvious that both elements influence detection of entanglement. Let us point out 
that the choice of maps is not limited to physical channels. The problem is usually 
formulated as follows: Given a density matrix of a bipartite system how strongly the 
two subsystems are entangled. That is we have a complete knowledge of the elements 
of the density matrix and we are allowed to execute arbitrary operation (function) on 
the matrix to calculate the entanglement. Such operation can be non-physical and even 
non-linear. Construction of entanglement witnesses using a general class of non-physical 
but linear maps has been investigated in Ref. p3]. The authors have showed that with 
the help of linear maps it is possible to distinguish perfectly the set of entangled states 
from the set of separable states. Here we show that this approach could be useful not 
only for the problem of detection but also for the problem of quantifying entanglement. 
Let us express the most general bipartite two-qubit state pu using the Pauli 
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operators aj, j = 1, 2, 3 

3 33" 

j=l k=l j;k=l 

where a^, (3k and 7^^ for j,k = 1,2,3 are real parameters. Further, consider a hnear 
map A12 

e 3 

A12 : P12 P12 + t(-1 + ^ifc^i ® ^k) ■ (21) 

j,k=i 

With the help of the map ( 12T1) and the trace distance we define the following function 

^(pi2) = -Tri2|Ai2[pi2] - P12I - 1 • (22) 
e 

The factor 1/e is there to eliminate the dependence on the epsilon while the —1 has 
been added for a convenience only. The function JF has the following properties: 

(i) ^(EjAjpia) < Ej AjJ'(pi2), convexity. 

(ii) (F{Ui®U2 P12 Ul®U2) = J-'ipu), local unitary equivalence Vt/i and Vf/2. 

(iii) J^ipu) > 0, Vpi2, non-negativity. 

(iv) J^i^pu) = 0, for all separable states. 

(v) (F{p\2) = C(pi2), where pi2 is pure or Werner state (for definition of the Werner 
state see Ref. [H]) and C(pi2) is the concurrence (see Ref. [15]). 

Through the extension of the proposed method to non-physical maps and a suitable 
choice of the map acting on the joint state pi2 we have managed to construct a function 
that detects entanglement on all Werner states. Moreover, some of the listed properties 
of the function JF are supposed to be fulfilled by a function that not only distinguishes 
separable and entangled states but performs a harder task - measures entanglement 
between two quantum systems. Though, the constructed function is not a proper 
measure of entanglement (there are entangled states for which JF is zero) a suitable 
extension might "correct" the function so that all entangled states are detected. 

4. channels 

In many physical situations it is less convenient to divide the system under consideration 
into two large subsystems than into a large number (say A^) of smaller (but equal 
systems). Typical example is the envisaged quantum computer composed of small 
micro-traps each holding a single qubit. In such case individual qubits are spatially 
separated so that the interaction with environment can be described by local maps 
Aj where the index i labels the qubits (or micro-traps). These maps can be derived 
phenomenologically or determined experimentally so their knowledge can be assumed. 
Obviously, we want to keep the infiuence of the environment as small as possible so each 
of the maps would satisfy a condition similar to Eq. ([T]) 

A(Ai[pi], Pi) < e ; = 1, ..AT , Vpi G S{Hi) , (23) 




Bounds on action of local quantum channels 



12 



where Di{., .)'s are again metrics (distance functions) and SiHi) is the set of all density 
operators for each i = 1, . . . N. Using these maps we can find out the state of a particular 
qubit after interaction with the environment. However, what is more important is the 
final state of the whole system 

Ai® Ajv[pi..7v], (24) 

and, in particular, how much the joint state pi„N has changed due to the interaction with 
the environment. This change can be characterized by a distance between the original 
state pi„N and the output of the product of the individual maps given by Eq. 

Di..;v(Ai ® . . . ® ^n[Pi..n], Pi..n) , (25) 

where the Di ^ is a metric (distance function) defined on the set of all density operators 
S{Ti-i,,N) of the joint system l..A^. Here we use the same definition of the metric 
(distance) as before and define the functions -Dj(.,.) for j = 1,...A^ and Di 
with the help of the norm of the corresponding Hilbert-Schmidt space (for more details 
see Sec. |2]) 

Di{p,a) =\\p- a\\i , (26) 
Di..N{p,or) = \\p- a\\i„N ■ (27) 
Using these definitions it can be shown that the distance in Eq. ( !25l) is always less or 



equal to A^y2 — 1/de where d is the dimension of the Hilbert space Hi. 

We note that the action of the product of local channels Ai . . . (g) Aat on separable 
states is such that Di. jv(Ai(8>. . .®Aiv[pi..Ar], Pl.at) < Ne. This means, that the restriction 
to the set of separable states leads to the decrease of the bound on Di by the factor 



To prove the statement we will use a very similar line of reasoning as in the case 
of two subsystems. First, taking advantage of the triangle inequality we bound the 
distance in Eq. as follows: 

^i..7v(Ai ■■An[pi..n],Pi..n) < 

Di..iv(Ai ® ..Aiv[pi..iv], 1 ® A2 ® ..An[pi..n]) 

+ Di„Nit ® ..1 ® An[Pi..n], Pi..n) ■ (28) 
Each of the terms on the right side of the last equation can be rewritten as 

Di„Nit® -t^Ai®!® .■t[p?,N],Pi!N) , (29) 

where 

p?^N = 1 ® ..1 ® Ai+i ® Ai+2 ® ..An[pi..n] , 

* Here we have used the bound ([7]) for separable states that can be easily extended to a multi-partite 
case. 
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so it is sufficient to bound the expression fl29l) . Next, we divide the whole system into 
two parts, an elementary system i and the rest. From this point the proof takes the 
same lines as in the case of two subsystems discussed in the Sec. El Therefore we recall 
the result Eq. (120|) obtained there and refer the reader to the Sec. |3] for more details. 
The Eq. (EO]) states that 

Di..7v(l® ..1® Ai®l® ..l[pPiv]>pPiv) < \/2- 1/c?, 

where d is the dimension of the i-th elementary subsystem. Since we have terms in 
the expression on the right in Eq. (1251) the distance fl2^ is bounded by 

D^„n{Ai ® . . . ® AnIpl.n], Pi..n) < Np-l/d e , (30) 

where A^ is the number of elementary subsystems each satisfying the condition (l26l) and 
d is the dimension of the Hilbert spaces Tii corresponding to the elementary subsystems. 

Example 

To illustrate the character of changes induced by the local maps on the global state of 
the whole system let us consider a simple model of A^ qubits undergoing a process of 
decoherence. That is the Hilbert spaces Tii are two-dimensional and the maps Aj are 
chosen to be 

K : ^(1 + a.a) ^\{^ + ^s^s + (1 - A;)K(Ti + ctsas]} , (31) 

where k is equal to k = \/2e in order to fulfill the conditions Eq. (l23l) . The action of the 
map Aj is such that it preserves the diagonal elements in basis formed by the eigenvectors 
of (Js while the non-diagonal elements are suppressed. Such maps describe the process of 
dephasing, a particular case of decoherence, since the vanishing of off-diagonal elements 
results into states that describe statistical mixtures. 

Consider the initial state of the joint system to be the Greenberger-Horn-Zeiliner 
(GHZ) state 

Pi..^ = ^{|0...0)(0...0| + |0...0)(l...l| 

+ |1...1)(0...0| + (32) 

The action of the map Ai ® ..A^v on the state Pl.at described above can be evaluated 
straightforwardly and we obtain 

Ai ® . . . ® An[pi..n] = ^ {|0 . . . 0)(0 . . . 0| + |1 . . . . . . 1| 

+ (l-A:)^(|0...0)(l...l| + /i.c.)} . (33) 

Despite the fact that the state of each individual qubit remains unchanged (a 
consequence of this is that the conditions (l23l) are trivially fulffiled) the state of the 
whole system changes because the off-diagonal elements are strongly suppressed. The 
distance ( l25l) between the input Pl.at and the corresponding output Ai ® . . . Kn[pi..n] 
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gives 

D(Ai ® .. ® An[pi..n],Pi..n) = [l + {l-kyN -2il-k)^] , 

which for e being very small can be estimated as 

D(Ai ® .. ® An[Pi..n], Pi..n) ~ ATe . (34) 

The deviation of the GHZ state under the action of the direct product of local maps 
Aj for sufficiently small e scales as A^e which confirms our more general result Eq. (l30ll . 
Though the result may seem to be optimistic (one might expect worse scaling with A^) 
the effect of the action of local maps is to disentangle the qubits (destroy quantum 
correlations between the qubits). In addition, the disentanglement itself is strong since 
the off-diagonal elements are suppressed exponentially with the increase of the number 
of systems involved in the dynamics. This example nicely illustrates that though the 
deviation expressed with the help of the distance fl25|) scales as A^e the entanglement 
may be destroyed much more dramatically. 

Finally note that in this example the bound for separable states A^e, derived with 
the help of Eq, ([7]), is not violated in spite of the fact that we have used an entangled 
state. We have already pointed out that it is not necessary for any entangled state 
to violate the separable bound. To show that the bound can be violated indeed one 
can choose the map Ai defined in Sec. |3] for maps Aj and the initial state of the form 
|6e//)®^/2 ^here \bell) denotes one of the Bell states (see, for instance. Sec. [3]). 



5. Conclusion 



We have analyzed the direct product of linear maps that describe local actions of a 
set of quantum channels. We have found a bound on the action of such product of 
maps (expressed as a distance between an input and output of the product) provided 
the linear maps composing the product are bounded as well. We have addressed two 
typical scenarios. In the first, a quantum system is divided into two subsystems and 
the product is composed of two maps acting on the two subsystems, respectively. In 
the second scenario a joint system is composed of A^ equal subsystems and we have A^ 
linear maps acting on A^ subsystems of a given quantum system. 

Our analysis has shown that the fundamental difference between the set of separable 
and entangled states yields two different bounds. For separable states the distance 

Eq. ([3]) is bounded by ^2 + 2^(1 — l/di){l — l/d2)e while in the case of entangled 



states the distance can be larger and is bounded from above by 2y2 — 1/dt. 

Let us note that the bound for separable states (fTTj) is optimal. That is there 
exists a pair of channels Ai and A2 such that the bound is saturated (examples are 
presented in Appendix C). It is interresting to note that the channels that saturate 
the bound are the same channels that saturate the separable bound ([7]) for the case of 



the trace distance (see Appendix C ) or the bound ( l20i) for entangles states in case of 
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two-dimensional systems. Clearly, to establish the upper bound on the action of a pair 
of local quantum channels it is sufficient to find a pair of channels for which the action is 
maximal and set the bound to this maximum. The form of the channels may depend on 
the dimensions di and d2- However, our results suggest that the channels for which the 
bounds are maximal are of the same form for arbitrary di and d2 and are the channels 
that we have used in our examples. 

In the end let us point out that our analysis is not restricted to the case of physical 
maps only and can be extended to the case of linear and hermicity preserving maps that 
are not physical. 
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Appendix A. 

We prove the relation 

||A[|A:)(/|]-|A:)(/| 1^ = ||A[ |/)(A:| ] - 1^ = 

\{mm{l\ + \l){k\)]-{\k){l\ + \l){k\)\\' 

+ ||A[(|fc)(/|-|/)(A:|)]-(|fc)(/|- 10(^1) ir}. (Al) 

for all physical (linear, hermicity preserving and completely positive) maps A with the 
norm defined in Eq.([2]) and k ^ I. 

Let us denote by V^i the expression A[|A;)(/| ] — |A;)(/|. Using the definition of the 
norm in Eq. ([2]) we have that for any physical map A 

\\A[\k){l\]-\k){l\ 
||A[ 10(^1 ]-|0(A:| IP 
l|A[(|A:)(/| + |/)(A:|)]-(|fc)(/| + 10(^1) 
\\A[{\k){l\-\l){k\)]-{\k){l\-\l){k\) 

Eq. (lA.ip is a direct consequence of the last result. 
Appendix B. 

In this appendix we prove two relations 

II A[(|fc)(/| + 10(^1) ]-(|A:)(/| + |/)(A:|) 
\\A[{\k){l\-\l){k\)]-{\k){l\-\l){k\) 



ti 



TT{VM + V^i)iV^i + V, 



kl 



kl 



< 2e, 

< 2e, 



(B.l) 
(B.2) 
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for all physical (linear, completely positive and hermicity preserving) maps A satisfying 
the condition given by Eq. ([1]) and k ^ I. The two expressions 

\\^[m{l\ + \l){k\)]-m{l\ + \l){k\)\\. 

\\H{\k){i\-\i){k\)]-{\k){i\-\i){k\)\\, 

can be rewritten as 

II A[ (pi -P2) ] - (Pi -P2) II , 

II A[i(p3 -P4) ] - «(P3 -P4) II , 

where 

Pi = \{\k) + mh.c.), 
P2 = \m-mh.c.), 
p, = \{\k)+tmh.c.), 
p, = \{\k)-imh.c.). 

By using the triangle inequality 

1 1 A[ (pi - P2) ] - (pi - P2) 1 1 < 1 1 A[ pi ] - pi 1 1 + 1 1 A[ p2 ] - p2 1 1 , 

1 1 A[ i(p3 - pi)]- iip3 - P4) 1 1 < 1 1 A[ P3 ] - P3 1 1 + 1 1 A[ p4 ] - P4 1 1 , 

we obtain the relations Eq. (1B.1I) and Eq. (IB. 21) owing to the conditions Eq. ([T]). 



Appendix C. 

Here we present an example showing that the bounds ([7j) and (fTTj) are optimal. In this 
example we will consider a more general case of distances Di{., .), D2{., .) and Di2{., .) 
and define the distances with p-norms 

DS(p„crJ = (Tr|p,-a,n^/P , 

where a labels the system 1, 2 or 12, p^ and cTq are density operators and p is a positive 
integer. The map Ai is chosen to be a contraction of the form 

Ai[pi] = (1- A;i) pi + ki -3-1, 

"1 

where di is the dimension of the Hilbert space Tii and ki is the contraction parameter. 
In what follows we assume fci = e/[(l - l/di)P + {di - l)/dlYlP so that the condition 
([T]) is fulfilled. Similarly, the map A2 is a contraction 

A2[p2] = (1-A;2) P2 + A;2 , 

"2 

where d2 is the dimension of the Hilbert space I-L2 and ki = e/[(l — l/(i2)^+((i2 — l)/'^!]^''^- 
For the initial state we choose a pure state of the form pi2 = |00)(00|. Keeping only 
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terms of the order of e the distance between the input pu and the output Ai (g) A2[pi2] 
gives 

D?2(Ai ® A2M,Pl2) = 



+ 



{di-l) + 



d2 



p \ i/p 

(4-1) . (C.l) 



Case study: Trace distance 

The trace distance is defined with the help of the 1-norm so that p = 1. The two 
contraction parameters ki and /c2 read ki = e/[2(l — l/(ii)] and /c2 = e/[2(l — i/d2)]- 
Using these relation in Eq. fIC.ip the distance between the input and the output reads 

Di2(Ai® Ai[pi2],Pi2) = 2e. (C.2) 



Case study: Hilbert- Schmidt distance 

The Hilbert-Schmidt distance is defined with the help of the 2-norm so that p = 2. 
The two contraction parameters ki and ^2 read ki = e/ ^Jl — 1/di and /c2 = e/ ^1 — 1/^2- 
Using these relations in Eq. (IC.ip the distance between the input and the corresponding 
output reads 



Di2(Ai® Ai[pi2],Pi2) = y2 + 2/(r^lM)(l-l/d2ye. (C.3) 
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